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' *\ ' We propose a Ramsey interferometry experiment using an entangled state of N atoms to reach the Heisenberg limit for the estimation 
r-y ' of an atomic phase shift if the atom number parity is perfectly determined. In a more realistic situation, due to statistical fluctuations of 
' the atom source and the finite detection efficiency, the parity is unknown. We then achieve about half the Heisenberg limit. The scheme 
involves an ensemble of circular Rydberg atoms which dispersively interact successively with two initially empty microwave cavities. The 
scheme does not require very high-Q cavities. An experimental realization with about ten entangled Rydberg atoms is achievable with 
, state of art apparatuses. 
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Quantum limits to noise in spectroscopy and interferometry jnH^EHEl have attracted an increasing 
£^ ■ interest in the last years. The precision of an interferometric phase measurement is ultimately limited by 
the "Heisenberg uncertainty relation" between energy and time p[J|BliniE3] ■ However, this limit can only 
Q_ ( ' be reached when N entangled particles are used. This represents a major experimental difficulty which 
has been overcome only very recently. In Ref. |Hj three beryllium ions, and more recently in Ref. (Hj six 
Ch \ beryllium ions have been prepared in a maximally entangled state and used in a Ramsey spectroscopy 
^ ■ experiment achieving a sensitivity close to the Heisenberg limit. In Refs. |11|I12| experiments involving 
q-i respectively three and four maximally entangled photons are reported, clearly showing interference fringes 
three and four times narrower than those obtained in an experiment with uncorrelated photons. Here, we 
propose an easy to implement cavity scheme, involving significantly more particles (N > 10). 
The aim of a typical interferometry experiment is the measurement of the relative phase <p between 
' two quantum states \g) and |e) of a two level atom. The standard approach is the Ramsey interferometry 
technique JB], in which a first tt/2 pulse produces a superposition of \g) and |e). After a waiting time, the 
relative phase is probed with a second n/2 pulse. The Ramsey setup is formally equivalent to an optical 
Mach-Zehnder interferometer, with the two tt/2 Ramsey pulses corresponding to the input and output 
50 — 50 beam splitters. The atomic phase shift corresponds to the differential phase shift between the two 
arms of the interferometer The uncertainty Aip of this phase measurement is limited by counting 

statistics to = l/y/N, where N is the number of detected particles (atoms for Ramsey spectroscopy, 
photons in the Mach-Zehnder interferometer) [7j. 

This limit holds, however, only if the N particles are uncorrelated, i. e. in a separable quantum state. 
If, instead, the N particles are quantum correlated and share some entanglement (which is equivalent to 
say that they are spin-squeezed |5lHl]L quantum noise sets a fundamental limit to the phase uncertainty, 
Atp = 1/N. This ultimate limit is called the Heisenberg limit and is achieved when a maximally entangled 
state of the N particles is used and appropriately measured [71ITU]. 

We present here a variant of a Ramsey interferometry experiment which approaches the ultimate 
Heisenberg-limited resolution and which is realizable with present technology. It involves a collective in- 
teraction of N Rydberg atoms with two high-Q microwave cavities. In the limit of a large atom-cavity 
detuning (dispersive interaction), the A^ atoms are prepared in a maximally entangled state, i.e. in an 
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atomic Schrodinger cat state after their collective interaction with the first cavity This state is sen- 
sitive to a Stark phase shift tp applied between the two cavities. The interaction with the second cavity 
reads out this phase shift. The final measurement of the atomic population reveals interference fringes iV 
times narrower than those achievable using non-entangled atoms. 

The outline of the paper is as follows. In Sec. [21 the scheme of the experiment is described in detail. 
In Sec. |31 we calculate the interference signal resulting from averaging over the fluctuations of the atom 
number. In Sec. |1] we present a modified scheme, which improves this averaged signal. Finally, in Sec. [5] 
a conditional detection scheme further improving the sensitivity of the detection of the phase shift is 
illustrated, while Sec. HO contains concluding remarks. 



2 Scheme of the interferometry experiment 



A schematic description of the proposed interferometry experiment is shown in Fig. ^ It involves a col- 
lection of N two-level atoms simultaneously crossing a set of microwave cavities. The two relevant atomic 
levels are circular Rubidium Rydberg states |e) and \g) with principal quantum numbers 51 and 50, and a 
transition frequency u eg /2ir = 51 GHz. The circular Rydberg states are characterized by a very long life- 
time (~ 30 ms) and a very large dipole moment |16| . The atoms are simultaneously prepared in zone B by 
the excitation of a velocity-selected atomic beam effusing from an oven. The N atoms sample, whose size 
is negligible at the scale of the millimeter wavelength A = 2irc/ui eg , crosses the arrangement of microwave 
cavities at the same velocity v. Three "Ramsey zones", Rj (j = 1,2,3), are separated by two identical 
microwave Fabry- Perot cavities C\ and C2, which have a high quality factor Q. The three Ramsey zones 
are microwave cavities with low Q, in which the Rydberg atoms can be subjected to classical resonant 
pulses, generated by standard phase- locked microwave sources. The three cavities Rj are resonant with 
the g <-> e transition, while the high-Q cavities C\ and C2 are off-resonant. The atom-cavity detunings are 
5i = uj eg — u C; i (u> Ci i is the mode angular frequency of cavity i). The atoms are finally detected in D by a 
state selective field-ionization detector, able to count both the number of atoms in |e), N e and in \g), N g . 

The ./V atoms collectively interact with the quantized radiation mode of C\ and C2, where they are 
symmetrically coupled The interaction of the atomic system with cavity i, in a frame rotating at co eg , 
is described by the Tavis-Cummings Hamiltonian JH] 

Hi,i{t) = -hSi&lai + ih^^- (at J_ - at J+) , (1) 

where we have introduced the angular momentum Dicke operators J_ = \g}i(e\, J+ = (J-)^ and 




Figure 1. Schematic description of the proposed cavity-QED interferometry experiment. The N atoms are prepared in the circular 
Rydberg state \g) in B, then cross the three classical field zones Rj and the two high-Q cavities Cj, and are finally detected in D. The 
squares denote it/ 2 Rabi pulses in Rj, the circle denotes the tunable Stark shift pulse in R2, and the diamonds denote the dispersive 

atom-cavity interaction in Cj. 



3 



J z = (1/2) Y2i \ e )i( e \ ~ lfi')j(fl , l I n Eq. 6j is the annihilation operator of the quasi-resonant mode 
in Cj, and £2j(t) = J7o ex P [— f 2 i 2 /u> 2 ] is the vacuum Rabi frequency describing the atom-field interaction. 
The time dependency is due to the Gaussian profile of the radiation mode with waist w by the moving 
atomic sample. 

We are interested in the dispersive limit in which the two cavities G\ and C2 are detuned far from the 
atomic transition, i.e., 5 3> QqVn, j c (jc is the cavity decay rate). Using cryogenic techniques together 
with appropriate "cooling" atoms ^B] > ° ne can approach the zero temperature limit and therefore we shall 
assume a mean thermal photon number of zero from now on. Neglecting the trivial constant of motion 
J 2 = (AT/2) (AT/2 + 1), Eq. (JTJ) is approximated as 

^disp,W = -^^(j z 2 -J 2 ), (2) 

explicitly showing that, in the dispersive limit, the atomic excitation J z and the photon excitation aja^, 
independently, are constants of motion. The term J 2 in Eq. © is responsible for the generation of atomic 
entanglement and spin squeezing [5]. This is reminiscent of the optical Kerr effect, which is characterized 
by a quadratic term in the photon number. In the atomic as well as in the optical case, this quadratic 
term is able, for particular values of the interaction time, to generate "Schrodinger cats", i.e., quantum 
superpositions of coherent states with different classical phases |15U19ll2*Uj . 

We adopt the compact representation of an atomic coherent state with orientation on a Bloch sphere 
(0, 0), (0 < < 7T, < (j)< 2vr), introduced in [2T], 




where the so-called Dicke states | J, m) are the eigenstates of J z , and J 2 with respective eigenvalues m and 
J(J+1). 

The atoms are initially prepared in the atomic coherent state |^)o = Y\i I d)i = 1 0, 0) and are then 
subjected to a sequence of unitary operations Ur. {i = 1,2,3) and Uc (j = 1,2) to generate the final 
state l^finai (see Figs. Hand 01): 

I lj>) final = Ur 3 Uc 2 Ur 2 (cp) U Cl U Rl I V) • (4) 

As a first step, the atoms undergo in Ri a it/2 pulse, described by the unitary operator = 
exp{iirJ y /2}. It aligns the collective Bloch vector along the x-axis, generating the atomic coherent state 

1^)^=^1^)0 = 1^/2,0). (5) 

The interaction with classical microwave pulses only changes the orientation of an atomic coherent state on 
the Bloch sphere |21j but the iV atoms remain disentangled. The atomic entanglement is produced in C\ by 
dispersive interaction with the radiation mode, which is initially prepared in the vacuum state. As shown 
in |T2], for specific values of the interaction time, an initial atomic coherent state \0,(ft) is transformed 
into a superposition of a finite number of coherent states with the same 9 and equally spaced values. 
An atomic Schrodinger cat state with two components is the maximally entangled A^-partite GHZ state 
corresponding to the superposition of two atomic coherent states with opposite orientations on the Bloch 
sphere, [|#, 0) + e^|7r — 9,<ft + n )] /y/2. In our proposal we obtain such an atomic cat state when ^5] 
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Figure 2. Schematic description of the state of the N Rydberg atoms on a Bloch sphere at the various stages of the interferometry 
experiment, both for even N and for odd AT. For both parities and at every stage, the N atoms are in a superposition of atomic 
coherent states (see Eqs. lsl- 1131 and Eqs. 1171 - 1181 for the explicit expression of the atomic states in the two cases). 



When this condition is satisfied, the interaction in C\ is represented by the unitary operator 



Ucy = exp i'~{Jl ~ Jz 



7T 



(7) 



and the state \ip)d = Uc 1 \iP)r 1 a t the exit of C\ becomes ^3] 
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if N is odd. This means that we have a cat state directed along the y-axis for even N and along the x-axis 
for odd N (see Fig.©. 

The goal of the interferometry experiment is to detect a variable phase difference (p between the two 
components of the cat state which is applied in R2 by a Stark pulse. This pulse induce a dephasing between 
the states \g) and |e), aligned along the z-axis. Thus, we have to rotate the cat state, aligned along the 
y-axis, to the z-axis prior to the application of the Stark pulse and then rotate it back into its initial 
direction. For the case of even N, this corresponds to the pulse sequence 



which is equivalent to Ur 2 (ip) = ex.p{iipj y }. We obtain 



e i7rJ x /2 e ii/?J z e -i7rJ x /2 



(10) 
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Thus, the two state components have acquired opposite phase shifts ±N<p/2, proportional to N. 
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Then the atoms enter C2, where they undergo the same transformation as in C\, Uc 2 = Uc 1 - It is easy 
to see that the atomic state becomes: 



\ even 
/C 2 



sin( 



N 

|5,7T> + (-1) 2 COS 



l|,o). 



(12) 



It is still an entangled atomic Schrodinger cat state, directed along the x-axis, with 99-dependent ampli- 
tudes. Since we always detect atomic energies, i.e., the number of atoms in \g) or |e), we use a final tt/2 
pulse in zone R3 before the detection, Ur 3 = Ur ± . The cat state of Eq. (fT2*|) is rotated around the y-axis 
and put along the z direction, so that one obtains the final state 



\ even 
/ final 



sin(^ 



0> + (-l)T cos(^) 



7T . 



(13) 



A measurement of the expectation value of (J z ) is done by counting the number of atoms in |e) = \9 = tt) 
and \g)= \0 = 0). Using the corresponding operators N e = N/2 + J z and N g = N/2 — J z , we obtain 



J z 



N e -N n 



(14) 



The expectation value of J z in state 
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— cos{Nip). 



(15) 



This signal displays Ramsey fringes with unit contrast oscillating N times faster than those corresponding 
to a single atom (see Fig. |3J). To calculate the sensitivity of (J z ) to changes of <p, we evaluate [22] 



\d(Jz)/d<p\ 



(16) 



where AJ Z = y { Jj 2 ) - {J z ) 2 measures the fluctuations of the operator J z . Eq. (O yields (J 2 2 ) even = N 2 /4 
so that, using Eqs. (|15j) - H16[) . one obtains Aip = 1/N, independent of ip. This means that, for an even 
number of atoms, we obtain a increase of sensitivity with respect to a standard Ramsey experiment 
involving N uncorrelated atoms |7j and thus achieve the Heisenberg limit. 

Grouping in a different way the unitary operations, we stress the link between this scheme and a genuine 
Ramsey interferometry experiment (see Fig. We divide the sequence of operations into three parts 
U{, (i = 1,2,3). The first, U\ = exp {— iirJ x /2}Uc 1 UR l , consists of the tt/2 pulse in R±, the dispersive 
interaction in C\ and the first tt/2 pulse in R2, (see curly brackets in Fig. These three operations 
prepare the N atoms, initially in the atomic coherent state \9 = 0), in the maximally entangled superpo- 
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'2 (in the case of N even). The second part corresponds to the 

application of the tunable phase (p through the Stark shift pulse, U2 = exp{i(/? J z }. Finally, the third part, 
^3 = Uji 3 Uc 2 exp {iirJ x /2}, consists of the second n/2 pulse in R2, the dispersive interaction in C2 and 
the final ir/2 pulse in R3. The operations U\ and C/3 are thus "Super-Ramsey pulses", corresponding to 
the tt/2 pulses of a classical Ramsey experiment. 

In the preceding discussion we only considered the case of an even atom number. For odd N values, the 
state \iP)r 2 at the exit of R2 becomes 
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for \(p\ < tt/2. For larger values of \<p\, similar expressions can be derived. 
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Figure 3. Interferometric signal ( j z ) after the entire Ramsey pulse sequence, according to Eqs. 1151 and 1191 , Shown are the cases for 
even (TV = 10, solid line) and odd (TV = 9, dashed line) atom number. The signals differ significantly in shape depending on the parity 

of the atom number. 

entangled cat state, however oriented in the xz plane at an angle (p with the x axis. The Stark shift pulse 
does not apply a relative phase on the state components, but rather rotates them around the z axis by an 
angle ip. As a consequence, the final state of the N atoms after R3 is very different from that of Eq. Q13|): 
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\(f,7r) + \<p, 0) - i(— 1) 2 Mtt — (p, 7r) — \tt — tp, 0} 



(18) 



It is a superposition of four coherent states, all in the xz plane (see Fig- 0) - Using the properties of atomic 
coherent states, it is possible to see that this state yields an interferometric signal given by 



(Jz 



\odd 



N , \ TV— 1 

— (COS if) 



(19) 



very different from that for even N (see Fig. EJ). It does not display any fast oscillation, but only 
dips at (p = 0, ±7r, having a width inversely proportional to y/N '. Moreover Eq. (fTH)) yields (J 2 2 ) odd = 
(iV 2 cos 2 ip + iVsin 2 y?) /4 so that one obtains the following phase sensitivity in the case of odd N 



Atp' 



odd 



yjN 2 [1 - (cos <p) 2N - 2 } - N(N - 1) sin 2 <p' 
N(N — 1)| sin ip\ I cos ^l^ -2 ' 



(20) 



This expression is minimum at tp = 0,±7r, where Aip odd = 1/yN, implying that, for an odd number of 
atoms, the phase sensitivity is always worse, or at best equal, to that of a standard Ramsey experiment 
involving iV uncorrelated atoms 



3 Interference signal averaged over the fluctuations of the atom number 



An experimental implementation of the scheme described above will be affected by the Poissonian fluc- 
tuations of the number of atoms and a non-ideal detection efficiency of the detectors, which makes it 
impossible to infer the parity of N. If we assume that we prepare iV Rydberg atoms with a Poissonian 
probability distribution with mean number N and average all detection events, the resulting interferomet- 
ric signal is given by the average over the distribution of a signal given by Eq. (|15|) for even N, and by 
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Figure 4. Interferometric signal {Jz) for Poisson distributed number of atoms, with mean number N = 10. The dashed lines show 
Eq. 1211 . valid in the dispersive limit. The full line is the numerical solution of the exact dynamics out of the dispersive regime, (a) 
8/Qq = 9.27, v — 40 m/s. (b) <5/f2o = 2.06, v = 180 m/s (see text for the other parameters). 



Eq. (jT9*)) for odd N. The corresponding expression is evaluated analytically: 



{Jz) = ye ^|sinh(ti 



^jcosf^Jcos^ 



cosh(n ¥ ,) [1 + sin(i^) sin ip] j , 



(21) 



with Up = Ncosip, Vu, = N sin <p. The behavior of this signal for N = 10 is shown in Fig. Eta) (dashed 
line), where it is compared with the prediction of a numerical solution of the dynamics driven by the 
exact Hamiltonian [Eq. in the presence of cavity damping (full line). The fast oscillating fringes of 
Eq. (|15j) except the central one at (p = are degraded by the incoherent average over different frequencies. 
In addition, the contribution of the experimental runs with an odd atom number reduces the height of the 
central fringe. 

The agreement of Eq. (|21|) with the numerical simulation is good since we have chosen a parameter 
region corresponding to the dispersive regime of large detunings. Using the values of Ref. |16j . i.e., S7o = 
0.31 x 10 6 s , w = 6 mm, an atomic decay rate 7 a = 33.3 s _1 and a cavity decay rate j c = 10 3 s , 
Fig-E^a) refers to = 9.27, i.e., 5 = 2.87 x 10 6 s _1 , which, taking the condition of Eq. © into account, 
implies an atomic velocity v = 40 m/s. The resulting signal {J z ) is periodic with period n. Its contrast is 
smaller than 1/2 because of the signal suppression in the wrong parity case. Despite that, this signal still 
scales with the number of atoms as the Heisenberg limit, thanks to the peak at <p = 0, whose width is 
inversely proportional to N, as it can be seen from the approximate expression of Eq. (|21j) at small (p 



{Jz) - -y exp I — ) sin 



<p(N + l) 



(22) 



The exact phase sensitivity is provided by Eq. (|16ft . whose explicit expression in the case of Poisson- 
distributed number of atoms can be obtained from Eq. (|21|) and the Poisson average of the squared 
momentum 



O Q AT O 

1 + cos (p + e sin cp 



(23) 



The resulting expression is a cumbersome function of cp, achieving its minimum at \ip\ ~ tt/(2N + 2), which 
also approximately corresponds to the points where the absolute value of the slope of the interferometric 
signal is maximum. At these phase shift values, A93 ~ 2/N, showing that at small values of (p and in the 
dispersive regime, our scheme reaches 1/2 of the Heisenberg limit despite the presence of a fluctuating 
number of atoms. 

The low atomic velocity is hardly compatible with the state of the art experiments |16j . We thus consider 
the case of a more realistic velocity v = 180 m/s, for which we expect to have a reasonable flux of atoms. 
In this case the 7r/2-pulse condition of Eq. © imposes a smaller detuning 5/VLq = 2.06. The behavior of 
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(J z ) in this non-dispersive regime is shown in Fig. 0(b) (full line), where it is compared with the dispersive 
limit expression of Eq. (|21|) for N = 10 (dashed line). The fringes at (p = ±7r have been washed out, 
while the central fringe at <p = is still visible, even though its contrast is now significantly reduced. The 
corresponding resolution in the estimation of <p is now decreased and one has Aip ~ 3.3/ N. 



4 Modified scheme with inversion of time evolution 

As shown in Fig. Efb), the contrast of the interference fringes degrades as soon as we leave the dispersive 
limit of slow velocities and we consider the fluctuations of the number of atoms. However, it is possible to 
significantly improve the scheme if we avoid the suppression of the signal at (p = in the odd N case (see 
Eq. (|19l0 which is responsible for the fact that the signal is practically never larger than zero. This can 
be achieved if the time evolution in C2 is "inverted", i.e., Hipif) — > —Hi^{t)- In the dispersive limit, this 
inversion can be achieved by changing only the sign of the detuning 82 through an appropriate Stark shift 
of the atomic levels. In such a case at (p = nothing happens in R2, and the two interactions in C\ and 
C2 cancel each other. The atoms only undergo the two ir/2 pulses in R± and i?3 which rotate the state 
from \9 = 0) to \9 = 7r). The final result is (J z ) = N/2 (N e = N, N g = 0) at <p = either in the even and 
in the odd case. Within the dispersive limit, we derive the state of the N atoms in the presence of this 
time evolution inversion. Nothing changes for even N, i.e., Eqs. (|12|) - (|15|) remain valid, while in the odd 
iV case the final state is now given by 



iodd,inv ' 

' final 2 



\<p,o) - - i(-iy 



TT 



ip,x) + |vr - ip,0) 



(24) 



yielding an interferometric signal 
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(25) 



Using Eq. (|25|) . the resulting Poisson-averaged interferometric signal can be obtained by changing only the 
sign of the term cosh (u^) in the expression of Eq. H21|). so that 



(Jz 



N 



-N 



sinh('u (/ ,) cos^) cos (p + cosh(n^) [1 — sm(vp) sin ip] 



(26) 




Figure 5. Interferometric signal ( j z ) mv in the presence of the "inversion" of time evolution in C2, for Poisson-distributed number of 
atoms with mean number TV = 10. The dashed line refers to Eq. 1261 . which has been derived assuming the dispersive Hamiltonian of 
Eq. while the full line refers to the numerical solution of the exact dynamics in the case of large detuning, S/Qo = 9.27 and 

v = 40 m/s (see text for the other parameters). 
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Figure 6. (a) Normalized interfcrometric signal (2J z ) lnv /N in the presence of the inversion in C2, for Poisson-distributcd number of 
atoms, with mean number N = 2, 4, 6, 8, 10, 12. The curves show the numerical solution of the exact dynamics in the non-dispersive 
regime, with 8/Q0 = 2.06 and v = 180 m/s. (b) Phase uncertainty Aip as defined in Eq. 1161 for the same parameter values, versus 
1/N. Stars refer to the minimum phase uncertainty and the dotted line is the least-squares linear fit of the data, while triangles refer to 
the phase uncertainty at the point of maximum absolute slope of the signal and the full line is the least-squares linear fit of the 

corresponding data. 

where again = N cos cp, = N sin tp. The behavior of this expression for N = 10 is shown in Fig. \5\ 
(dashed line), where it is compared with the corresponding value of {J z ) mv obtained from the numerical 
solution of the dynamics driven by the exact Hamiltonian of Eq. and in the presence of cavity damping 
(full line) for the same parameters chosen in Fig. ^a). 

The signal is again periodic with period equal to tt, but it now shows clear peaks at <p> = 0, ±7T, where 
it reaches the maximum value N/2. The width of the peaks is again inversely proportional to N, because 
Eq. can be approximated at small ip as 



a> inv ^yexp 



Nip 2 



cos 



(p(N + l) 



(27) 



The explicit expression of the phase sensitivity of Eq. (|16|) can be obtained from Eq. 1)26|) and from Eq. (|23j) 
for ( which is valid also in the presence of the inversion. The resulting expression is again a cumbersome 
function of <p, now achieving its minimum at \ip\ ~ ir /(AN + 4), where Aip ~ 1.45/iV. Differently from the 
preceding scheme without inversion, the point of minimum phase uncertainty does not coincide with the 
point of maximum slope of the signal (J z ) mv , which is again approximately equal to \<p\ ~ ir/(2N + 2), 
as in the scheme with no inversion of time evolution of the preceding Section. At these values of ip, the 
phase uncertainty is larger, but still scales as the Heisenberg limit with the number of atoms, because it 
is Aip ~ 1.65/iV. 

The inversion of the time evolution in C2 becomes of crucial importance in the non-dispersive regime. 
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In this case, one also has to change the sign of J± in addition to reversing the sign of 62 in order to get 
Hi 2,(£) — > —Hi 2(t). This second inversion can be realized by applying two opposite tt Stark shifts. A first 
one, described by the operator exp{i7r J z }, must be applied in R2 just before the entrance in C2, while the 
second, described by the operator exp{— iirJ z }, must be applied in R3 soon after the exit from Ci- The 
new unitary operators for the Ramsey zones R2 and R3 then read: 

U'^) = ^U R M: (28) 
U' R3 =U R3 e-™ 3 *. (29) 

As a result of the inversion, the narrow central interference peak at <p = becomes much more robust and 
it survives even outside the dispersive regime of slow velocities. This behavior of (J z ) mv is shown in Fig. El 
for different values of N, while in Fig. E^b) the dependency Acp ~ 2.0/ N is clearly visible. This fact proves 
that thanks to inversion in C2, we are able to reach about one half of the ultimate, Heisenberg-limited, 
resolution for the estimate of ip, even in the non-dispersive regime of fast atoms. 

Notice that in the non-dispersive regime the dynamics is no more described by the Hamiltonian of Eq. (J2J) 
and therefore the atomic phase shift is no more given by Eq. Nonetheless, we have still used Eq. © 
for the choice of the value of 5 at a given velocity v, because it turned out in the numerical simulation to 
give the best results. When fast atoms are used (v = 180 m/s), the atomic state after the interactions in 
C\ and C2 is rather different from that of the dispersive limit [Eqs. (jSJ and l|12|)]. However, the interference 
experiment in the presence of the inversion in C2 still approaches the Heisenberg limit. 

The main experimental limitations affecting this scheme are cavity damping, the spontaneous emission 
of the atoms and the non-unit detection efficiency Cavity damping is not a serious limitation because the 
cavities play a passive role in the experiment, only allowing the virtual exchange of excitations between the 
atoms. The cavities are initially in the vacuum state and can get photons from the excited Rydberg atoms 
only close to resonance. Cavity decay starts to significantly affect the generated maximally entangled GHZ 
atomic state only when j c ~ 5 ~ flo^/N. This is not a strong condition on j c because typically 5 ~ 10 6 
s _1 . Thus, a moderate quality factor of Q = u> e g/jc > 10 6 would be sufficient, compared to experimentally 
attainable values of Q > 10 8 . 

Also atomic spontaneous emission does not represent a serious problem, thanks to the use of circular 
Rydberg levels with lifetimes of the order of 30 ms, even though it puts a practical upper bound on the 
maximum number of atoms N. Spontaneous emission implies a decoherence timescale of the maximally 
entangled atomic GHZ states of the order of 30/iV ms. This decoherence time has to be longer than the 
time of flight through the apparatus which is of the order of 1 ms. Therefore it is reasonable to perform 
the experiment with up to 20 atoms. 

To consider a non-ideal detection efficiency, we assume that the detector is characterized by the same 
quantum efficiency 77 for detecting atoms in e or g, that there are no dark counts, and we neglect the 
possibility of a wrong state detection (e *-+ g), the phase uncertainty simply reads A<// = r/ -1 Ac/2. Typical 
detection efficiencies attainable in a cavity QED experiment with Rydberg atoms are between 80% and 
100% H3] . 



5 Conditioning to the number of detected atoms 

In the preceding schemes, we did not use any information on the number of detected atoms. We now 
show that we can improve the phase sensitivity by means of an appropriate data processing of the signals 
conditioned to the number of detected atoms. If one post-selects only the data corresponding to a given 
number of detected atoms = N^ et + Ng Bt , in the ideal case of a detector with quantum efficiency rj = 1, 
one would get an interference signal coinciding with the ideal one [Eq. (|15|)] with N = (even), thanks to 
the conditional generation of a maximally entangled GHZ state of atoms. If instead the detector is not 
perfect, the interference signal is the result of an average over many experimental runs in which the actual 
number of atoms is equal or larger than N^. As a consequence, the contrast of the interference fringes 
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rapidly worsens for a decreasing detection efficiency r\. This post-selection strategy is in a certain sense 
opposite to the one considered in the preceding sections, where we have averaged the interferometric signal 
over all the runs, regardless the value of the number of detected atoms. In the post-selection case, one 
gives zero weight to the signal conditioned to a number of detected atoms Nj, different from the selected 
one. The Poissonian-averaged signal of the preceding section instead corresponds to keep all the data and 
to give the same weight to the conditioned signals. 

It is evident that the post-selection strategy is not optimal for achieving the maximum phase sensi- 
tivity, because it implies wasting most of the resources, i.e., all the atoms of the runs with a number of 
detected atoms different from the selected one N d . In this experiment, the best strategy for achieving the 
maximum phase sensitivity is to keep all the data, however grouping them into different sets according 
to the corresponding number of detected atoms A^, and then give appropriate weights w(Nd) to the sig- 
nal conditioned to the detection of atoms, ( J z (Nd)) cond ■ This corresponds to consider the following 
interferometric signal 



(J z ) w = P(N d )w(N d )(J z (N d )) 
N d =o 



cond 



(30) 



where P(N d ) is the probability of detecting N d atoms which, due to the assumptions made above, is given 
by 



P(N d ) = e 



-T)N 



(31) 



The optimized signal corresponds to take the weights w opt (N d ) which, for each <p, minimize the phase 
uncertainty Aip of Eq. (|16ft , which for the signal of Eq. (|30)1 has the explicit form 



Aip 





i) 2 (J*(N d )) cond ~ (T,Z=o p (N d )w(N d )(J z (N d ))^ d f 


1/2 




E% d ^P(Nd)w(N d )d(J g (N d ))°<»«/dip 





(32) 



The Poisson-averaged signal of the preceding sections corresponds to the particular case w(N d ) = 1 VN d 
in Eq. (j3U|) . 

We illustrate this optimization strategy by applying it to the experiment with the inversion of time 
evolution in C2, in the dispersive limit of large detuning. In this limit, using Eqs. (|15|) and (|25|) . it is 
possible to find, after long but straightforward calculations, that the conditional signal is given by 



(J z (Nd)) C ° nd = { e R > cosip [cos(7V> + N x simp) + (cosc^ 1 ] 



+(-l) nd e 



N d -Ni cos tp 



[ca&(N d (p - Ni simp) - (cos^)^ -1 ] j 



(33) 



where iVj = A^(l — 77), is the mean number of atoms lost by the detector. Moreover in the dispersive limit 
one also finds 



(4 2 (^)> C ° nd 



A^ 2 N d -N 2 d . 2 



+ 



sin cp 



Na-N, 



-l) JVd e 



(34) 



These two latter expressions are then inserted into the expression for the phase uncertainty Aip of Eq. 1)32(1 , 
which is then minimized with respect to the weights w(N d ). The corresponding optimal weights w opt (N d ) 
depend upon the phase shift applied in R2, tp, due to the dependence of Aip upon </?, and when inserted 
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-7t -7t/2 tc/2 7C 5 10 15 20 M 25 

cp (rad.) iN d 



Figure 7. (a) Optimal interferometric signal ( j z ) opt of Eq. 1351 (full line) versus the phase shift applied in R2, ip. The optimal 
weights are determined by the condition of minimum phase uncertainty at each (p. The dashed line refers to the Poisson-averaged signal 

of Eq. 1261 . corresponding to take equal weights, (b) Weights w opt (N l i) of the optimal signal versus at the point of the maximum 
phase sensitivity. The plots refer to rj = 0.8 and TV = 12.5 and to the dispersive limit of large detuning, where Eqs. 1331 and 1341 apply. 

into Eq. Q3UJI . one gets the optimal interference signal 

oo 

(Jz)° pt = E P(N d )™° pt (N d )(l(N d )) COnd - (35) 

N d =0 

This optimized signal is shown in Fig. Efa) (full line), which refers to the case r] = 0.8 and N = 12.5. The 
corresponding minimum phase uncertainty is a function of (p, and its minimum value A(p m [ n is closer to the 
ideal Heisenberg limit: Ay?min — 1.3/r]N. Comparing with the result of the preceding section, we see that 
optimizing the information provided by conditioning on the number of detected atoms yields an improved 
phase sensitivity, since the corresponding minimum phase uncertainty of the Poisson-averaged case scales 
as A9J m i n ~ 1A5/t]N . As it happens in the Poisson-averaged signal, the best phase sensitivity is achieved 
at a phase shift (f smaller than that corresponding to the maximum slope of the optimal signal {J z ) , 
which is approximately equal to \(p\ ~ ir/(3r]N). At this value of the phase shift, the phase uncertainty 
still scales as the Heisenberg limit, but is larger and it is roughly given by Atp ~ 1.5/r/N. 

It is also interesting to see the behavior of the optimal weights w opt (N ( i) at the phase shift <p corresponding 
to the maximum phase sensitivity, which are shown in Fig. Efb), again for the case r\ = 0.8 and N = 12.5. 
One can notice that they monotonically increase for decreasing N d and this is due to the fact that the 
optimal weights tends to be inversely proportional to ( (N d )) cond . From the analysis of Section II, one 
could have expected an oscillating behavior of w opt (N d ), with maxima corresponding to N d even and 
minima corresponding to N d odd. However it can be seen that this happens only in the limit ry — * 1, when 
the probability to detect a wrong parity of the number of atoms becomes negligible. Only in this limit, 
one can safely discriminate between an even and an odd number of atoms and suppress the contribution 
of the runs with an odd number of atoms. 



6 Conclusions 

We have presented a Ramsey-like interference experiment for a cavity-QED system, able to reach the 
ultimate Heisenberg limit for the estimation of an atomic phase-shift. We considered a system of N two- 
level Rydberg atoms successively crossing two microwave cavities. A dispersive atom cavity interaction 
is able to generate an atomic Schrodinger cat, i.e., a superposition of two atomic coherent states, which 
represents a maximally entangled state of the ./V atoms. Using this state, we designed an interference 
experiment yielding fringes N times narrower than those one would have obtained if N disentangled 
atoms were used. 

As discussed in Sec. |3 the most important limitations affecting the the experiment are the fluctua- 
tions of the number of Rydberg atoms, and the non-unit detection efficiency Despite these limitations, 
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one can approach Heisenb erg-limited sensitivity, because a narrow central fringe, with a width inversely 
proportional to N, survives even when averaged over the fluctuations of the number of detected atoms. 

Finally we have considered a conditional scheme in which we post-select only the events with a fixed 
number of detected atoms. In such a case, the best strategy for achieving the maximum phase sensitivity 
is to keep all the data, however grouping them into different sets according to the corresponding number 
of detected atoms Nd, and then give appropriate weights w(Nd) to the signal conditioned to the detection 
of Nd atoms, {J z {Nd)) cond ■ One can then determine the optimal weights maximizing the phase sensitivity 
and the corresponding optimized signal closely approaches the ultimate Heisenberg limit, because the 
minimum phase uncertainty is A<^ m j n ~ 1.3 /rjN. 

The proposed experiment should allow for the first time to achieve the Heisenberg limit for spectroscopy 
with a larger number of atoms compared to the experiments realized so far with entangled photons or ions. 
Despite this application in spectroscopy, characterizing the maximally entangled Schrodinger cat state by 
e.g. state tomography and monitoring its decoherence would be a very interesting perspective. 
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